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' Abstract 
Oh 

^ We discuss the interplay of low-x physics and QCD scahng violations by extending the unified 

approach describing inclusive structure functions and diffractive production in j*p interactions 
proposed in previous papers, to large values of Q^. We describe the procedure of extracting, from 
the non-perturbative model, initial conditions for the QCD evolution that respect unitarity. As- 
^ ^ suming Regge factorization of the diffractive structure function, a similar procedure is proposed 

for the calculation of hard diffraction. The results are in good agreement with experimental 
data on the proton structure function F2 and the most recent data on the reduced diffractive 
cross section, xpar ■ Predictions for both F2 and Fl are presented in a wide kinematical 
y—( range and compared to calculations within high-energy QCD. 

o 
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1 Introduction 



The interaction of a highly energetic virtual photon with a nucleon, or a nucleus, probes the high- 
energy limit of strong interactions. In this situation it is important both to find the correct degrees 
^ of freedom of the interaction and to preserve unitarity of the cross section. In the infinite momentum 

frame the latter comes about, in the limit of large gluon density, as non-linear gluon interactions 
which tame the growth of parton distributions and eventually leads to saturation. In the target 
rest frame the same phenomenon is manifested in terms of multiple scattering. 

Experimental data indicate a fast growth of the 'j*p cross section accompanied by a large 
diffractive cross section at the highest available energies [T-21 . This calls for a unified treatment 



of several effects: a proper unitarization of the total cross section at low momentum scales in 
conjunction with a correct treatment of QCD scahng violations. While present day data give a rich 
insight into the latter, the question of whether saturation effects have been observed is still under 
debate. All the same, these effects are of crucial importance when extrapolating to higher energies, 
deep inelastic scattering (DIS) off nuclei and heavy- ion collisions. 
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DIS measurements are typically used to constrain universal parton distribution functions (PDFs) 
through an analysis of QCD scaling violations. The most recent state-of-the-art fits based on the 



DGLAP evolution equations up to NNLO have been presented in 122 , 23 . Regrettably, one has no 



perturbative theoretical control of the initial conditions and, thus, extrapolations to kinematical 
regions yet unexplored by experiments are affected by significant uncertainties, e.g. as obtained by 



recent PDF fits using neural networks 24 



An established approach to the problem of high-energy scattering in QCD is the dipole model 
|25[|26| . One can show that at very high energies the 7*p scattering can be modeled as a convolution 
of the wave function of a virtual photon fluctuating into a qq dipole and the subsequent dipole- 
proton cross section. While the former is known exactly from QED, the latter can be treated 
theoretically (approximated e.g. by two-gluon exchange) or simply parameterized. Still, the dipole 
model contains several weaknesses which limit its range of applicability, such as lack of scaling 
violations owing to QCD evolution. Furthermore, the model is only valid for small dipole sizes and 
at low-x. 

Using perturbative arguments one can show that the strength of the qq — N cross section should 
scale with the transverse area of the dipole 25|26| . A particular simple parameterization of this cross 
section was sugg ested in [27]. The so-called GBW dipole cross section encodes the characteristic 
features of parton saturation by assuming an x-dependent damping of large-size dipoles, or, in 
other words, a saturation scale oc x~^. A fit to experimental data resulted in a value A ~ 0.28. 



The model was extended to describe diffraction in 28 . 

In fact, the effective energy dependence of the cross section, Ae// = dlnF2/dln{l/x), changes 
from a quite small value at low-Q^, consistent with the so-called soft pomeron, to a steady growth 
in the perturbative region. The dipole model incorporates some scaling violations in the 7* — qq 
wave function which mimic the QCD behavior, but fails to match with the DGLAP equations at 



high-Q^ data due to the lack of scale dependence of the qq — N cross section. The authors in 29 
attempt to include these by making use of the connection between the dipole cross section and 
the gluon distribution function at leading logarithmic accuracy. At the initial scale they include 
solely a non-zero gluon PDF and perform a LO QCD fit to constrain its parameters. An impact 



parameter dependence was also introduced in a later extension of the model 30 . Calculations of 
diffraction within the improved dipole model have also been presented 

Recently, a general evolution equation for high-energy QCD has been derived using renormaliza- 



tion group techniques, the so-called JIMWLK equation 32-38 



accuracy, by the significantly less complicated BK equation 38 , 39 . Unfortunately, calculations at 
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leading order result, for realistic values of the strong coupling constant, in an x-dependence of 
the saturation scale significantly larger than the A extracted phenomenologically |27|. Sub- leading 
effects, such as energy-momentum conservation and inclusion of running of the strong coupling 



constant, are believed to tame this growth 40,41 . It was not until recently that a complete nu- 



merical solution 42 of the BK equation with running coupling was presented [43j-[45j . In [46] a set 
of parameters related to the initial distribution was fixed to obtain a good agreement with F2 and 
Fl at X < 0.01 and predictions down to x = lO"^'^ were made for a wide range of Q^. 

In spite of the successes, there are many aspects of the theoretical computations which are 
still to be fully understood, such as the impact parameter dependence and the description of 
diffraction. In the absence of a unified QCD approach to the entirety of 7*A^ (7*^) processes, a 
useful guidance for investigating the connection between non-perturbative and perturbative aspects 
of DIS valid for extrapolations to extremely small momentum fractions, x, can be found in terms of 
the Reggeon calculus f47| with a supercritical pomeron, Ajp = ajp — 1 > 0, and the partonic picture 
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of 7*A^ interactions. Multi-reggeon exchanges are included so as to satisfy s-channel unitarity. In a 
particular realization of these models, in the 7* wave- function, one distinguishes explicitly between 



a large (L) and a small (5) component 48 , 49 . The former interacts strongly even at high-Q'^ but 
is quite rare, while the latter interacts according to r ~ l/Q- Thus both components have a leading 
dependence of the total cross section while the L-component gives the leading contribution to 
diffraction (1/Q^ vs. 1/Q^). Large-mass diffraction is included through triple-reggeon interactions. 
In this work we will follow the approach of |49|, where the small component was cast in the form 
of the dipole model. The model gave a simultaneous description of inclusive F2 and diffraction in 
the region of < ^ 5 — 10 GeV^, and was used to predict the structure functions at very low 
X. In jSOj these results were used to predict nuclear shadowing calculated in the Glauber-Gribov 
theory, in good agreement with data. In Sec.[2]we give a detailed description of the model (formulae 
for diffraction are given in Appendix [A|) . 

With the advent of high-energy colliders, such as HERA and LHC, and the planned electron- 
hadron experiments 5T||52 , the need for low-x structure functions for nucleons and nuclei at high-Q^ 



have arisen. This motivates an extension of the model mentioned above [49] to the perturbative 
regime by the inclusion of QCD scaling violations. We describe a prescription for extracting the 
initial conditions at leading order for the DGLAP equations from the non-perturbative model 
both for inclusive F2 and diffraction. In the former case, this procedure does not involve new 
parameters. The situation for the inclusive diffractive cross section is more complex, because it 
involves both more complicated reggeon exchanges and additional variables in the problem. For the 
proper description of data in the whole (3 and xjp region we identify explicitly pomeron and reggeon 
contributions to diffraction. One can then invoke a supplementary factorization of variables, the 
so-called Regge factorization |53|, which allows for a comprehensible QCD analysis. In the reggeon 
case, we identify and include missing diagrams which are crucial for a proper description of data. 
Details on the initial conditions and subsequent QCD evolution in the inclusive and diffractive cases 
are given in Sees. [3] and [5} respectively. 

Thus, equipped with properly unitarized initial conditions for the DGLAP evolution equations 
we obtain leading-order structure functions and PDFs for the proton down to x ~ 10~^ at high-Q^. 
The resulting F2 and xpFj^j^ are shown to be in good agreement with the most recent experimental 
data. We also compute the longitudinal structure function within the dipole model using the 
perturbative gluon PDF thus obtained in Sec.[4| Comparisons are made with the recently computed 
solution of the running-coupling BK equation |46j. Finally, we present our conclusions in Sec.[6j 

2 Brief description of the CFSK model 

At small-x, the total cross section for 7*p interactions is related to the structure function F2{x, Q^) 

by 

cx?;*)(H^^Q^) = ^^^F2(x,g2), (1) 

where x = Q"^ / {W'^ + Q^), W = ^/s is the invariant mass of the produced hadronic system and 
is the virtuality of the photon. In the model of [49], denoted in the following as CFSK, the total 
cross section was written as a sum of two contributions 

= j dHa^^f^{h^s,Q'), (2) 
4? (6, Q') = gUQ>t\b, s, Q') + at\b, s, Q') . (3) 



3 




The first term on the right hand side of eq. ([3]) describes the non-perturbative interaction of large- 
size partonic configurations of the virtual photon with the target, while the second term describes 
the interaction of small-size configurations. The function g'j^{Q'^), determining the coupling of the 
7* to the large qq dipole, was chosen in the form 



gi(o) 
1 + qyml 



diQ') = rf^^ , (4) 



vanishing at high while the Q^-dependence of the small-size component is an inherent charac- 
teristic of the dipole model. 

The cross section of the L-component in impact parameter space is cast in the quasi-eikonal 
form (to include proton dissociation) 

ar(6,.,Q^) = ^---H-CXL{b,s,Q^)) 

where the function XL accounts for multiple pomeron (IP) and reggeon (IR = {/, A2,... }) exchanges 
and triple-reggeon interactions, as follows 



1 + ax3(s,b, Q^j 



(6) 



The first term in eq. ^ corresponds to a summation of fan-type diagrams, also called the Schwim- 
mer model 54 . The constant a is defined as a = 5r^(0)rjpjpjp/167r, where 5^(0) is the proton- 
pomeron coupling and rjpjpp is the triple-pomeron coupling. The eikonal functions Xlo ~ IP,IR) 
are written in the standard Regge form 



where 

A, = a,(0)-1 , e = ln^^^ , X(i, = Rl,L + c^'k^ (8) 

•So + Q 

and ak{0), a'^ are the intercept and slope of the corresponding trajectories, respectively. The 
function ^ is chosen such that it behaves as ln(l/x) at high-Q^ and ln(s/so) at = 0. The 
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triple-reggeon interaction term, X3is,b,Q'^) is 



1 



Ajp (xjp) 

I3max (j^p 



exp 



4 Ajp {x]p) 
X3{s,b,Q\f3) , 



\X]pJ 



(9) 
(10) 



where, in this case, X]p{x]p) = R^p + In (l/j;jp), where xp = x//? is the Bjorken momentum 
fraction of the pomeron, and the limits of integration are given by Pmin = xjxp^f^^ = lOx and 
fimax = Q'^/iQ'^ + ^Imin)- triple-reggcon eikonal accounts for heavy-mass diffraction, thus 

Mmin = 1 GeV A generic diagram of the 7*p-interaction is shown in Fig. [ij 

The cross section of the S'-component has been cast in the standard dipole form 



(tot) 



E 

T,L 



dr 



dz 



V'^'^(r,z)"a5(r,s,Q^), 



(11) 



where rg is a cut-off parameter on the size of the dipoles, to be fitted, and ■0^^^^ are the wave 
functions of the qq pair corresponding to transverse and longitudinal polarizations of the virtual 
photon, the corresponding squares given by 



ip {r,z) 



6ap 



47r2 



E 
E 



^2 + (1 - zfe'^Kfier) + m^^Kl{er) 
\q^z^ + (1 - zfKlier) 



(12) 
(13) 



Here = z{\ — z)Q'^ + rrig, and Kq and Ki are McDonald functions. We use the same quark 
mass, ruq, for all three quark flavors. The dipole-nucleon cross section as{r, s,Q'^) can be written 
in terms of the cross section at a given impact parameter 



^s{r,s,Q^) = 4j d%as{r,b,s,Q^) , 



(14) 



where crsir, 6, s, Q'^) is cast analogously to eqs. ^ and Q for (Tl(6, s, Q'^) with xl replaced by xs- 
The dependence of the dipole cross section on r is introduced taking into account that for small 
dipoles cross sections are proportional to r^, so that the eikonal function is defined as 



xfo(&,e) 



exp Ajp^ 



4A§p(e), 



(15) 



where = R^pg + a'jp£,. Note that secondary (iR) exchanges do not contribute to the S- 
component. 

A notable feature of the model is its growing interaction radius with energy, encoded in the 
functions Aq^i. and A^, which leads to an increase of the cross section at low x. It is also worth 
noticing that the amount of damping of the functions xl and X5) see eq. Q, are controlled by 

^In the CFSK model X3 contains a contribution from IfUPlP diagrams (exchanges in the triple-reggeon or 
pomeron diagrams are labeled clockwise starting from the left). These are unphysical and have removed in the 
calculations presented in this paper. Numerically, this contribution is insignificant in calculations of inclusive DIS 
and are only sizeable in the high-a;ip, low-/3 region of diffraction which we discuss in detail below. 
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Fixed parameters 


CFSK f49 


CFSK' 


lS]p 


n 9 
u.z 




4.56 X 10" 


-3 


O.DO X iU 




n ^ 

-u.o 




1.87 GeV- 


-2 


z.yo V 


a'jp 


0.25 GeV-2 


Cf 


0.56 GeV" 


-2 


0.46 GeV-2 


a'f 


0.9 GeV-2 


So 


0.79 GeV 


2 


0.79 GeV2 


R2 

^OkL 


3 GeV-2 


a 


4.63 X 10-2 GeV-2 


6.13 X 10-2 GeV-2 


r2 


2 GeV-2 


ml 


0.59 GeV 


2 


0.70 GeV2 


^Ik 


2.2 GeV-2 


Cs 


0.185 




0.105 


If 


8 


ro 


1.06 GeV" 


-1 


1.33 GeV-i 


c 


1.5 


ml 


0.15 GeV 


2 


1 X 10-3 GeV^ 



^ Put to -0.5 in the CFSK' fit. 
Table 1: Parameters of the 7*A^ model: the column labeled CFSK corresponds to the original 



model parameters, as found in 49 ; the column labeled CFSK' corresponds to a new fit performed 



with a smaller value of the reggeon intercept, Aj. 



the parameter a and the function X3 which are strongly constrained by diffraction. In this way, 
unitarization of the cross section and the amount of diffraction are intimately linked. 

Several of the model parameters are fixed from studies of the energy dependence of total and 
diffractive cross sections of hadronic reactions in Regge theory, while the rest were fitted to inclusive 



and diffractive DIS at low x and < < 5 GeV^ 49 . A summary of the parameter values is 
given m Table [ij the fitted ones from [49| can be found in the column labeled CFSK. This particular 
realization of the model conjectured a (non-standard) value of the reggeon intercept, Aj = —0.3. 
We have redone the fit to the same dataset fixing Aj = —0.5, which in the following we will denote 
as CFSK', see rightmost column in Table [T] for details. The new fit works even slightly better than 
the old one and has several additional advantages which we will describe in detail below. It is 
interesting to notice that the fitted cut-off on the size of the dipoles in the S'-component is about 
0.2^0.25 fm. 



3 QCD evolution of the inclusive structure function 

In order to generalize the CFSK model to large values of , QCD scaling violations have to be 
included. The model is therefore used as the initial condition for the DGLAP evolution equations 
at the initial scale = 2 GeV^ (any > 1 GeV^ can, in principle, be set as the initial scale 
and the results are rather insensitive to the value of Qo)- The initial condition for the evolution 
equations should be given for all x. While the CFSK model is valid only for small values of x, it 
can easily be extended to the x ~ 1 region. In order to do so, we follow a standard procedure [55] 
multiplying the partonic distributions by the relevant powers of (1 — x) as explained below. In 
what follows we will work at leading order (LO). 

Let us first infer the valence quark PDFs from the total cross section. They correspond to the 
exchange of secondary reggeons. In our model these are only inlcuded in the L-component, see 
eq. ([5]). In order to separate the iR-contribution, only linear terms in xff will be taken since it is 
sizeable only at not too low x, where multi-reggeon exchanges can be neglected. We therefore find 
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that the valence quark contribution to F2 from the model is 



r^low — X 

2V 



For the proton, we take n'^^'-^Cx, Q^)/2 = d''^'"-''{x,Q'^) [55] and use 



ilow—x I 



:plow—x 
2V 



9 



xu 



V 



9 



(16) 



(17) 



as the low-x valence quark distribution at the initial scale. The extension to high-x for the proton 
is carried out by multiplying the uy-quarks by (1 — x)"*^'^ and the dy-quarks by (1 — x)"^*^ 
respectively, where 



n(Q2 



3 
2 



1 + 



and c = 3.55 GeV 55 . Finally, the valence quark distributions are given by 



xny(x, Qo) 
xdv{x, QI) 



2R 



low—x , 



f: 



2V 

ilow—x , 



2V 



x,Ql) (1 - X, 
x,Ql) (1-x) 



n{Ql)+l 



(18) 



(19) 
(20) 



The sea quark PDF is given by the sum of the S-component and the singlet contribution to the 
L-component, i.e. neglecting all iR-terms in the latter: 



rploW — X 

2 Sea 



{tot) 



+ a 



(tot) 



c;=o 



(21) 



In order to obtain the PDF's for the different flavors in the sea, we define ^(x, Q^) = u = u = d 
d = 2s = 2s0so that 



11 



xS' 



low—x 



(22) 



Taking into account the relevant (1 — x) factor for the high-x behavior, the sea quark distribution 
is finally given by 



x5(x,Q^) = _^F^o--(^,Qg)(l_^ 



11 



(23) 



at the initial scale. 

At sufficiently large the term xs can be related to the distribution of gluons in the proton. 
In the leading logarithmic approximation 56 , 57 we have that 



(24) 



where xg{x, Q^) is the gluon distribution function of a proton and 05 is the strong coupling constant. 
For small values of xs^ we obtain 



xg 



low—x 



6 



cf 



dH 



1 + aX3(&,'S,Q2) 



(25) 



^Here we neglect the difTerence between u and d quarks in the region x < 0.1. 
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Figure 2: Initial quark and gluon parton distribution functions in CFSK (left) and CFSK' (right) 
models compared to the CTEQ6 LO parameterization [23] at Qq = 2 GeV^. 



The extrapolation to large values of x is the same as for the sea quarks with an additional (1 — x) ^ 
factor [55], that is 

xg{x, Ql) = xg^'^-'^ix, QI) (1 - x)"(Qo)+3 . (26) 



Thus, the equations (19)-(20), (23) and (26) constitute the initial conditions for the DGLAP evo- 
lution equations for the inclusive structure function. 

The treatment of heavy quarks is done in the zero-mass variable flavor number scheme. This 
scheme provides matching prescriptions between Nf and Nf + 1 evolved PDFs at a given threshold 
scale, fiT, proportional to the heavy quark mass, mg. The proportionality constant is not known 
theoretically but can be estimated requiring smoothness of observables, see [58[|59| . We take the 
charm quark mass to be rric = 1.4 GeV, neglect bottom and top contributions and take ht = 2.5 nic- 
Further investigation of the impact of heavy quarks and a detailed comparison to heavy quark 
structure functions, F2 and F2, lie beyond the scope of the present work. 

The parton distribution functions at the initial scale at leading order are thus fixed unambigu- 
ously as described above and they have to fulfill the valence and momentum sum rules, namely 

/ dx uv{x, Qq) = 2 , 

Jdx dv{x,Ql) = 1 , (27) 
/ dxx{uv{x,Ql) + dv{x,Ql) + hS{x,Ql) + g{x,Ql)^ = 1 . 

The partonic decomposition of the original CFSK model does not automatically fulfill these rules, 
which would force us to introduce uncomfortably large overall normalization factors. The CFSK' fit 
described above, on the other hand, automatically fulfills the sum rules to a good approximation. 



We have verified that eqs. (27) are satisified with an accuracy better that 5 10%. To further 
improve the sum rules we have increased the parameter C^^ by 12%, see Table [l] The fraction of 
the proton momentum carried by the gluons at the initial scale is then 

/ dxxg{x,Ql) _ f 0.59 in CFSK , 



Jdxx{S{x,Ql)+g{x,Ql)) [ 0.49 in CFSK' 



(28) 
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The initial parton distribution functions (PDFs) for both fits are compared to the CTEQ6 LO 
parameterization [23| in Fig. [2j The choice of Aj = —0.5 is clearly improving the CFSK' valence 
PDF contribution, and the smaller value of the Cg parameter than in the original CFSK brings 
the gluons closer to the fit of CTEQ6. Since the QCD sum rules seem to be better satisfied in the 
CFSK' model with no additional normalization of the input PDFs, we will continue using this set 
of parameters in what follows. 

We have subsequently performed QCD evolution at LO, with as{M'^) = 0.126 using the QCD- 
NUM evolution code [60]. The F2, and the associated parton distribution functions, thus obtained 
are denoted CFSK'e (evolved). Comparing with data from experiments [l|-[T6], a total of 847 data- 
points, the resulting x^/d.o.f. is 2 both for CFSK' (0 < < 2 GeV^) and in the evolved case 
(2 < Q2 < 1000 GeV2)0 The results of CFSK' and CFSK'e are shown and compared to a set of 
the experimental data in Fig. [3| Summarizing, we have shown that the proposed model is in good 
agreement with experimental data in the whole range of accessible and down to very low x. 

The LO DGLAP evolution takes into account large logarithms in Q^, yet at very high energies 
one expects large corrections arising from the missing terms in ln(l/x) which are more properly 
accounted for in the BK equation or in linear resummation schemes [61[|62| . In Fig. |4]we compare 
predictions of low-x F2 to the recent solution of the BK equation presented in |46| , where an initial 
qqN cross section in the form of the GBW model was assumed. We notice large differences between 
the models at very large Q^, which may partly be caused by the fact that the calculation presented 
in |46j does not take into account heavy flavors which are important at those momentum scales 
and partly by the absence of impact parameter dependence in the calculation |46| . Regardless of 
this, the deviations between the models for x < 10~^ at > 10 GeV^ hint of a breakdown of the 
DGLAP equations. Although the CFSK model satisfies unitarity, the observed difference is purely 
an effect of high-energy QCD and will be extremely interesting to study in the future. 



4 The longitudinal structure function 

The measurement of the longitudinal structure function of the proton at HERA has been eagerly 
awaited. The longitudinal structure function is zero at leading order QCD and is a direct measure 
of the size of the gluon distribution which had before this only been accessed through the scaling 
violations of the total structure function. Thus, it is believed that it can be a sensitive probe to 
saturation effects at low x. In the dipole model at leading logarithmic accuracy Fl only probes the 
dynamics of small, perturbative configurations. 

Although the longitudinal structure function is a NLO observable in the QCD improved parton 
model within collinear factorization, one can calculate this quantity in the framework of the dipole 
model using the connection between the dipole-nucleon cross section and the gluon distribution 
function at LLA. On the one hand, we can use directly the form of the dipole-nucleon cross section 



of the unevolved CFSK' model, eq. (|26|) and on the other hand we can make use of the relation 



of this cross section to the gluon distribution at > Qq, eq. (24), which we extract from the 



perturbative calculation, CFSK'e. 

Note that in the latter case there is a slight mismatch between the non-perturbative and pertur- 



bative gluons due to the integration over impact parameter. In principle, the definition in eq. (26) 
gives us the first term in a series of the (impact parameter dependent) quasi-eikonal model. To be 
able to resum the quasi-eikonal series we should include the correct impact parameter dependence 



^Here, we have not considered the normahzation errors of the data sets when calculating the 
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Figure 3: The structure function F2 in the CFSK'e model compared to data. The blue curve is 
the unevolved CFSK' model and the red solid curves is the CFSK' evolved at LO. Black points are 
experimental data (T-16 . 



of each individual term. We therefore define the impact parameter dependent gluon density in full 



10 




analogy with each of the pomeron terms as follows 



xg{b,x,Q'^ 



e^p{-b^/4Xljp{0} 



xg{x,Q 



where we have ensured the proper normalization. The final qq — N cross section is therefore 



(29) 



as{r,x,Q'^) = A / cfb 



2C 



1 - exp <^ -C 



xg{x, Q^) 



, (30) 



where both the gluon distribution and the strong coupling constant are calculated at LO. The 



resulting longitudinal structure function is found by convoluting as in eq. (30) with the wave- 



function for longitudinally polarized photons in eq. (13). 



We compare the two prescriptions for the dipole-nucleon cross section described above to HI [17| 
and ZEUS ^18j data in Figs. [5] and [6j At the present moment the experimental data have too 
large errors to distinguish between models, but more precise data at low-Q^ could increase the 
discriminating power of the observable. 

Finally, Fl at low x is expected to be sensitive to saturation effects which could be probed at 
a future electron-proton collider. The authors of [63j found that combined data on F2 and Fl give 
a strong discriminating power in revealing saturation effects. In particular, we note an order of 
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Figure 5: The longitudinal structure function of the proton as a function of in CFSK' and 
CFSK'e compared to data averaged over x from HI [Tt] (left figure) and ZEUS (right figure). 
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Figure 6: The x-dependence of the longitudinal structure function of the proton for several bins 



in the CFSK' and CFSK'e models compared to HI data 17 . 



magnitude difference between the prediction of the CFSK model and the BK model at low x and 
Q^, which are compared in Fig. [Tj This region lies beyond the reach of perturbation theory and 
can shed light on the transition to the non-perturbative regime. 
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Figure 7: Predictions for Fl in the CFSK model and from a numerical solution of the BK equation 
|46j. Note that the lower red, dashed curve of the BK solution corresponds to = 10~^ GeV^. 



5 QCD evolution of the difFractive structure function 



Using the AGK cutting rules 64 we can readily obtain the diffractive cross section from the formulas 
of the total cross section described above. Once more, the diffractive cross section will consist of 
contributions from the large and small partonic configurations of the virtual photon wave-function, 
while the leading contribution arise only from the former, unlike in the inclusive case. Therefore, 
the CFSK model includes also the explicit contribution from SIP diagrams which are responsible 
for high-mass diffraction (the low-/3 region). For specific details on the diffractive part of the CFSK 
model we refer the reader to [48 , 49 and Appendix [A} 



The collinear factorization theorem of the diffractive cross section, X]pa^^'^\f3,xip,Q'^,t), is 
valid at fixed xjp and t only for the resolved photon [65| . Nonetheless, experimental data show to 
a good approximation that diffractive DIS data satisfy proton vertex factorization, whereby the 
dependences on variables which describe the scattered proton (xjp, t) factorize from those describing 
the hard partonic interaction (Q^,/3). This property is also known as Regge factorization [53] . 
For example, the slope parameter B, extracted by fitting the t distribution to the form da/dt oc 



e'^*, shows no significant variations from the average value 19 . Also, one observes no significant 
variation of the pomeron intercept with 19 . 



These observations hold as long as one-pomeron exchange dominates the cross section. The 
CFSK model involves more complicated diagrams and respects this naive factorization only in a 
limited kinematical region. On one hand, the multi-reggeon exchanges change the effective intercept 
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Figure 8: Regge factorization in the CFSK model. 



of the pomeron. On the other hand, the low-/3 and high-Xip region is dominated by reggeon 
exchange, a missing piece in the original formulation of CFSK. For the sake of simplicity, comparing 
our calculations to the experimentally available t-integrated (reduced) diffractive cross section, 
X]pa^^^\p, xjp, Q^), we will make use of the general Regge factorization, representing the diffractive 
structure functions as 

f!§{xjp, /3, Q2) = -fjpi^{xjp) Ff (/3, Q2) + n^/K/pCxj.) Ff (/3, Q^) . (31) 

where F^ are the reggeon/pomeron structure functions and /^^^ is the t-integrated {k = IRi,IP) 
flux factor of the proton 

-I tent Xjp''^ ' 

where tmin = iT^j^jXp and tcut = — 1 GeV^. Finally, njj is an unknown normalization of the reggeon 
contribution. 

(3) 

Concentrating on the pomeron contribution to Fg-p for the time being, we have checked numeri- 
cally that the CFSK model of diffraction follows a universal trend, such that xjpFj^ x Xjp ~ f{l^) 
for a large range of xjp at the inital scale Qg, where ^eff = 1 — ciijp^ {^) = 0.123 is the effective 
pomeron slope, see Fig. [8] for details. Thus, Regge factorization is very well satisfied in the CFSK 
model, a slight breakdown observed at low /3 values and xjp > 10~^. 

The pure multi-pomeron contribution at a given xjpq can consequently be defined as 



:jpF^^ {xip, /3, g2) = Ajp xjpQ Ff Q^) , (33) 



where A]p is an overall normalization and the HUPIP contribution is not included. We will hereafter 
choose X]PQ = 0.01 as a reference. Thus, the value of the diffractive structure function at a given 
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xjp is simply 



xipF^jl{x]p,l3,Q'^ 



xiPo \ 



2A 



(3) 



2V 



(34) 



We carry on with an attempt of a partonic decomposition of the pomeron structure function, 
analogously to the inclusive case. 

The normalization of the pomeron flux factor, Ap, is unknown and therefore we define the 
singlet quark diffractive parton distribution at the initial scale as 



/35^(/3,Q^) 



9 
11 



xjpF^^{x]p,l3,Ql) 



(35) 



where we have assumed three active quarks flavors and suppressed the strange quarks by a factor 
2 in analogy to the proton PDFs explained in Sec. [Sj The brackets [...]q denote that this value is 
taken at a given xjpq . The general xjp-dependence is given through the relation in eq. ( 34 ) . 



Concerning the gluons, since we are not able to separate the flux and the diffractive PDFs 
(dPDFs), we effectively obtain the /3 and dependences through the QCD evolution and gain 
access to a product of the Pomeron flux and gluon dPDF given by 



X]pfiP/p{x]p) (3g^{f3,Q^ 



f X]Pq\ 
\X]p ) 



2A 



(36) 



where the normalization is the same for gluons and quarks. There is no a priori procedure to 
extract the gluon distribution from the CFSK model in the diffractive case. Thus, the last factor 



on the right hand side of eq. ( 36 ) is parameterized as 



^fiP. Q')l , = AgP^'^il - /3)^«(1 + Dg^) exp 



0.001 

1 — x 



(37) 



so that it is integrable. Since, in the single-pomeron exchange model, the ratio of gluon and sea- 
quark distributions in the proton should equal the corresponding ratio in the pomeron at low x 



and /3, respectively, we fix the parameter Bg to the standard value, Bg 



-Ajp. The remaining 



parameters, Ag, Cg and Dg, of the gluon dPDF at the initial scale have to be extracted from the 
data. 



The description of diffractive data (Q^ > 2 GeV^) 
a careful analysis of the necessary components 



in the whole region of xp and (3 demands 



66 . As shown in 20 , 671 , the inclusion of reggeon 



terms is crucial for describing data at high xjp. A missing piece in the original model is the HUPIR 
contribution which is dominant in the triple-reggeon region, i.e. at low /3 and high xp. Its diffractive 
cut corresponds to a reggeon-exchange with large-mass diffraction. One should accordingly intro- 
duce parton densities in the reggeon and, following the standard procedure 21 , we identify them 



with the pion ones; together with the standard expression for the reggeon flux, see eq. (32), we thus 



(3) 

obtain the reggeon contribution to i2X' • present work we have employed the LO DGLAP fit 

from l68j, which was shown to be in reasonable agreement with recent leading neutron data from 
HERAr|69|. The parameters of the reggeon trajectory and t-slope are the experimentally extracted 
ones, listed in 21 . Finally, the normalization of the flux was set so that xp //™™ fiR/pdt = 1 at 

Xp = XpQ. 

We have performed a LO DGLAP fit (with as{M'^) and fiT as described in Sec. [s]) of the 
gluon dPDF parameters and np of the reggeon contribution together with the input sea dPDF 
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A, 








fit result 


0.108 


-1.82 


-0.91 


0.0107 



Table 2: Parameters of the gluon density of the pomeron and the normalization of the reggeon 



contribution obtained from a fit to diffractive data 19 ,21 



given in eq. (35) to experimental data on the diffractive structure function from HI |2l] (LRG, 461 
points) and ZEUS 19 (LRG, 277 points and LPS, 118 points). The resulting x^/d.o.f. for both the 



unevolved and evolved CFSK' models with tt-jr = are large. The largest deviations arise at high 
X]p where precise data dictate a more careful treatment of the reggeon contribution to diffraction. 
Leaving njR as a free parameter we obtain a x^/d.o.f. of 1.8, improving the agreement significantly. 
The resulting values of the parameters obtained from the fit are listed in Table [2| 

The QCD evolved CFSK' initial condition with the fitted parameterization of the gluon dPDF 
and the reggeon contribution are compared to the HI LRG data for the lowest bins in Fig. [9] 
together with the unevolved CFSK' model. The overall description of the data seems satisfactory 
in both models, except the low-/? bins where the missing reggeon contributions leads to larger 
deviations for the model without evolution. The reason for the fairly good description of data at 
high in the original model can, in fact, be traced back to the dipole form of the small component 
of the 3iP contribution which gives rise to a logarithmic growth in . 

Whereas the description of the inclusive proton F2 at > 2 GeV^ calls for a extension 
of the CFSK model while the scaling violations of the diffractive cross section are much better 
accounted for, we have presented, for the sake of consistency, a complete procedure for a partonic 
decomposition of the original model and subsequent QCD evolution for both the inclusive and 
diffractive calculations, dubbed CFSK'e. The proper magnitude of the scaling violations of the 
diffractive part is crucial for the comparison with the inclusive F2 in the calculation of nuclear 
shadowing in Glauber-Gribov theory |70j, giving rise to the correct Q^-evolution of shadowing. 

Finally, as an additional cross-check of the region of validity of the CFSK'e model presented in 
detail above we study the condition that adi ff I otot < 1/2 related to the conservation of unitarity, 
also called the Pumplin bound (tTI. We define 



\dxFB{xp)F'ii{fi,XF,Q'') .... 
= • ^''^ 



In Fig. [To] we plot Rp for several 0.5 GeV^ < < 20 GeV^. We see that Rp is below 0.5 down 
to X ~ 10-12^ which proves the validity of the model for large in the kinematical region we are 
discussing in this paper, relevant for present and future experiments and even cosmic-ray physics. 

6 Conclusions 



We have extended the original CFSK model 48 , 49 of DIS to the high-Q region by the inclusion 



of QCD scaling violations via the DGLAP evolution equations. The extended model (CFSK'e) was 
used to calculate F2, Fi^ and xwF'^{x]p^ /?, Q^) of the proton for < < 1600 GeV^ and all values 
of X down to 10"^. The agreement with existing experimental data on all of these observables is 
good. 

The CFSK'e model for scattering off nucleons can be used for robust predictions at high energies 



for a wide range of observables, from DIS observables to multiplicities in heavy-ion collisions 72 
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0.5 




Figure 10: The Pumplin ratio for the CFSK and CFSK'e (for > 2 GeV^) models. The calculation 
for = 3.5 GeV^ is also shown for a larger range in x in the sub-view of the plot. 
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It also provides a bridge between perturbative and non-perturbative regimes of QCD and can serve 
as a baseline for deviations from the standard Q^'^volution involving resummations of logarithms 
of X. Finally, it bridges unitarity corrections with the strength of diffraction at high energies. 

We observed differences between our approach and high-energy QCD evolution encoded in the 
solution of the BK equation in kinematical regions beyond current experimental reach. At low x, 
F2 at high is significantly larger in our calculations. This may imply the breakdown of the 
collinear factorization leading to the DGLAP evolution and a transition to a linear resummation or 
non-linear regime. On the other hand, the BK calculation is performed at a fixed impact parameter 
while the CFSK explicitly includes a growing interaction radius with energy, which make it hard to 
quantify these deviations. At low Q^, Fi may shed light on the transition between the perturbative 
and non-perturbative (Regge) regimes. 

As we have pointed out above, the model for low-x structure functions at high respects 
unitarity at the initial scale of QCD evolution, but lacks nevertheless logarithms of 1 /x that should 
be resummed at high energies. One should also keep in mind that the QCD analysis presented above 
was done at leading order in the coupling constant. In the future we will attempt to include higher- 
order corrections in our calculations. An extension to the nuclear case is also under developement 
|73l. 
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A The diffractive cross section in the CFSK model 



The diffractive cross section of a virtual photon in the CFSK model 49 is given by three terms 



a. 



IP 



Eidiff) , 
i=L,S 



(39) 



where 



a 



(diff) 
L 

(diff) 
S 



(TJPIPIP 



(tot). 



/ / d 



r I dz 
Jo 



V'^'^(z,r) [as{r,b,s,Q') 

2\^-2CxL{b.s.Q^) 



d^b d\ dz i;^^'^{z,r) XripAb,s,Q^)e 
Jo Jo 



2\-2CxL(b.s.Q') 



(40) 

(41) 
(42) 
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Here 



xjpiP]p{b,s,Q'^) = a^{b,s,Q^)x3ib,s,Q'^ 



(43) 



where i = L,S, given by the first term in eq. ^ and X3 is defined in eq. (10). To find the 

(reduced) diffr active structure function as a function of (3, defined by 



X]pF^^{x]p,P,Q'^ 



dt xjp 



da 
dxpdt ' 



we perform the same decomposition as above, giving 



7(3) _ p(3) , p,(3) p(3) 



(44) 



(45) 



where the individual pieces are extended to the whole /3-region. In particular, the S-component 
amounts to 



where ^ = (Q2 + sq)/{Q'^ + M^) = px/x, and 



/3„ 



dl3_ 



(46) 



(47) 



with f3„ 



10 X and f3„ 



Q2/(Q2 _|_ 4,7j,2)_ Next, the L-component is 



^ ^(3) _QW)^ 



(0) 



dH (xf) AA[/3^»+^^-^^(1-/3) 



n(Q2) 



(48) 



c=o 



49 



where xf is given by the first term in eq. ([6]). Note the difference between our formula and 
where also double IPIR and IRIR exchanges were taken into account (numerically, these contributions 
are insignificant). Finally, the 3iP-component is 



■^P^2ViIP 



(0) 

J-jp-r2X)3jp 



47r2 



Of 



2a / d26x3(6,^,Q^/3) [4*°*)(s,6,Q2)+4*°*)'-^''^*(,,5,q2) 



(49) 



,(50) 
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